Theta vacuum and entanglement interaction 
in the three-flavor Polyakov-loop extended Nambu-Jona-Lasinio model 
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We investigate theta- vacuum effects on the QCD phase diagram for the realistic 2-1-1 flavor system, using 
the three-flavor Polyakov-extended Nambu-Jona-Lasinio (PNJL) model and the entanglement PNJL model as 
an extension of the PNJL model. The theta-vacuum effects make the chiral transition sharper. For large theta- 
vacuum angle the chiral transition becomes first order even if the quark number chemical potential is zero, when 
the entanglement coupling between the chiral condensate and the Polyakov loop is taken into account. We 
finally propose a way of circumventing the sign problem on lattice QCD with finite theta. 

PACS numbers: 11.30.Rd, 12.40.-y 



I. INTRODUCTION 

The existence of instanton solution in quantum chromody- 
namics (QCD) requires the QCD Lagrangian with the theta 
vacuum: 

Cqcd = qf{l^D,^ + mf)qf + iK'^F!lu 
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in EucHdean spacetime, where F^^^ is the field strength of 
gluon [1]. Theoretically the angle 9 can take any arbitrary 
value between — tt and tt. However, experiments indicate 
\9\ < 3 X 10"^" 10. HI- The Lagrangian is invariant under 
the combination of the parity transformation and the transfor- 
mation 9 — !■ —9, so that parity (P) and charge-parity sym- 
metry (CP) become exact only at 9 — and ±7r; note that 
9 = — TT is identical with 9 = tt because of the periodicity 27r 
in 9. Why should 9 be so small? This long-standing puzzle 
is called the strong CP problem; see Ref. (0] and references 
therein for the detail. 

For T higher than the QCD scale Aqcd, there is a pos- 
sibility that 9 is effectively varied to finite values depending 
on spacetime coordinates {t,x), since sphalerons are so acti- 
vated as to jump over the potential barrier between the differ- 
ent degenerate ground states |5]. If this happens, P and CP 
symmetries can be violated locally in high-energy heavy-ion 
collisions or the earl^ Universe at T « Aqcd- Actually, it is 
argued in Refs. [6, 7] that 9 may be of order 1 at the epoch 
of the QCD phase transition in the early Universe, whereas it 
vanishes at the present epoch lf8l- [l2ll . This finite value of 9 
could be a new source of large CP violation in the early Uni- 
verse and may be a crucial missing element for solving the 
puzzle of baryogenesis. 
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In the early stage of heavy-ion collision, the magnetic field 
is formed, while the effective 9{t, x) deviates the total num- 
ber of particles plus antiparticles with right-handed helicity 
from those with left-handed helicity. As a consequence of this 
fact, an electromagnetic cuiTent is generated along the mag- 
netic field, since particles with right-handed helicity move op- 
posite to antiparticles with right-handed helicity. This is the 
so-called chiral magnetic effect The chiral mag- 

netic effect may explain the charge separations observed in the 
recent STAR results |16]. The thermal system with nonzero 9 
is thus quite interesting. 

For vacuum with no temperature (T) and no quark-number 
chemical potential (/i), parity P is preserved when 6* = [ ITJ, 
but is spontaneously broken when 9 = n ifisi [l9ll . The P 
violation, called the Dashen mechanism, is essentially non- 
perturbative, but the first-principle lattice QCD (LQCD) is 
not applicable for the case of finite 9 because of the sign 
problem. Temperature (T) and/or quark-number chemical po- 
tential (/i) dependence of the mechanism has then been ana- 
lyzed by effective models such as the chiral perturbation the- 
ory |i7.,2Q-,24,1, the Nambu-Jona-Lasinio (NJL) model |,25.-28ll 
and the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) 
model f23,^. 

Using the two-flavor NJL model ifsil - O . Fujihara, Inagaki 
and Kimura made a pioneering work on the P violation at 
9 = n f23\ and Boer and Boomsma studied this issue exten- 
sively |26, 27]. In the previous works |29, 30], we extended 
the formalism to the two-flavor PNJL and entanglement PNJL 
(EPNJL) models and investigated effects of the theta vacuum 
on the QCD phase diagram. Very recently similar analyses 
were made for the realistic case of 2h-1 flavors by using the 
NJL model [28]. It is then highly expected that the finite-6' 
effect is investigated in the 2h-1 flavor case by using the PNJL 
and EPNJL models that are more reliable than the NJL model. 

If QCD with finite 9 is analyzed directly with LQCD, one 
can get conclusive results on the theta vacuum. Here we con- 
sider a way of minimizing the sign problem in LQCD. For this 
purpose we transform the quark field q to the new one q' by 
the following SUa{5) ® C^a(1) transformation. 
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The QCD Lagrangian is then rewritten into 



cqcd = ^i^M + ^s^^^s + 4^;.^;. (3) 

l—u^d 

with the new quark field g', where 

Mi[9) = -fuD^+micoB{e/2) + mii-i^s\n{e/2), (4) 
X, EE -f^D^+rUs. (5) 

Only the Dirac operator Aii{9) has 6' dependence in (|3]l. The 
determinant of satisfies 



AciMi{0) = (detA^K-^*))' 



(6) 



The sign problem is thus induced by the 9-oAd (P-odd) term, 
m;i75 sin (0/2). The difficulty of the sign problem is ex- 
pected to be minimized in the QCD Lagrangian of Q, since 
the 6'-odd term includes the light quark mass mi that is much 
smaller than /Iqcd as a typical scale of QCD. This point is 
discussed in this paper 

In this paper, we analyze effects of the theta vacuum on 
the QCD phase diagram for the realistic case of 2+1 flavors 
by using the three-flavor PNJL and EPNJL models fsl. 
Particularly, the three-flavor phase diagram is investigated in 
the ji-T plane with some values of 9. Through the analysis, 
we finally propose a way of circumventing the sign problem 
on LQCD calculations with finite 9. 

This paper is organized as follows. In Sec. |II] we recapit- 
ulate the three-flavor PNJL and EPNJL models. In Sec. |III1 
numerical results are shown. Section |IV] is devoted to sum- 
mary. 



II. FORMALISM 

The three-flavor PNJL Lagrangian with the 6'-dependent 
anomaly term is obtained in Euclidean spacetime by 

C =q{j^D^ + rha - ■^ifi)q 

8 

a=0 

+ Gd [e'*' det - 75)?^ + e'*'' det q,{l + 75)9^- 

L li li 



-U{<P[A],<P*[A],T), 



(7) 



where = ~ iSuiAIXa/l with the Gell-Mann matri- 
ces Aq. The corresponding PNJL Lagrangian in Minkowski 
spacetime is shown in Refs. ll29i[30ll . The three-flavor quark 
fields q = {qu, qd, qs) have masses rho = diag(m„, m^, m^), 
and the chemical potential matrix /t is defined by /i = 
diag(/x, fi, /i) with the quark-number chemical potential /i. 
Parameters Gs and Gd denote coupling constants of the 
scalar-type four-quark and the Kobayashi-Maskawa-'t Hooft 
(KMT) determinant interaction ll36i ISTIi . respectively, where 
the determinant runs in the flavor space. The KMT determi- 
nant interaction breaks the C/a(1) symmetry explicitly. Obvi- 
ously, the theta-vacuum parameter 9 has a periodicity of 2tt. 
We then restrict 9 in a period < 9 < 2tt. 



The gauge field Af^ is treated as a homogeneous and static 
background field in the PNJL model IM H HI HI HH^H. 
The Polyakov-loop and its conjugate 'P* are determined in 
the Euclidean space by 



'?=itrc(L), <?* = itrc(Z), 



(8) 



where L — exp{iA4/T) with A4/T — diag((/)r, 0^) in 
the Polyakov-gauge; note that Aq is traceless and hence (j)r + 
<l>q + (f'b = 0. Therefore we obtain 



3 
3 
3 



(9) 



We use the Polyakov potential U of Ref . 114111 : 



-<p*<p 



b{T) ln(l - 6'P'P* + A{<P^ + -P*^) ~ 3{'l><P*f) 



(10) 



with 



aiT)^ao+a,{^)+a,(^)\ biT) ^ b,(^y . 

(11) 

The parameter set in U is fitted to LQCD data at finite T in the 
pure gauge limit. The parameters except To are summarized 
in Table U The Polyakov potential yields a first-order decon- 
finement phase transition at T = To in the pure gauge theory. 
The original value of Tq is 270 MeV determined from the pure 
gauge LQCD data, but the PNJL model with this value of Tq 
yields a larger value of the pseudocritical temperature Tc of 
the deconfinement transition at zero chemical potential than 
Tc « 160 MeV predicted by full LQCD jsstissll. Therefore 
we rescale Tq to 195 (150) MeV so that the PNJL (EPNJL) 
model can reproduce Tc = 160 MeV iQ. 



ao 


ai 


12 


h 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential 
sector determined in Ref. f4r]. All parameters are dimensionless. 



Now the quark field q is transformed into the new one q' by 
dU in order to remove 9 dependence of the determinant inter- 
action. As shown later, this transformation provides the ther- 
modynamic potential i7 with a compact form and furthermore 
convenient to discuss the sign problem in LQCD. The present 
three-flavor PNJL model has 18 scalar and pseudoscalar con- 
densates of quark-antiquark pair, but flavor off-diagonal con- 
densates vanish for the s^tem with flavor symmetric chemi- 
cal potentials only ||26[ l27i l29i [soll . Since the quark-number 
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chemical potential considered in this paper is flavor diagonal, 
we can concentrate our discussion on flavor-diagonal conden- 
sates. Under the transformation (|2]i, the flavor-diagonal quark- 
antiquark condensates, c/ = q/q/ and 77/ = QfilaQf, are 
transformed into ct', 



q'^q'j: and 77^ 



Qfilsq'f as 



m 

OS 

Vs 



cos 
— sm 

V's 



(f)^;+sin(f)r/;, 



(12) 

(13) 
(14) 
(15) 



for I — u,d. The Lagrangian is then rewritten into 

-GsY^iiq'Xaqr + iq'^l^XaqT] 



a=0 

Gd [det qlil - -f5)qj + det q[{l + 75)9; 



with 



mo+ 



nio- 



diag(m„+,md+,ms+) 

diag (cos (f)m„, cos (|)m<j, m, 

diag(m„_,md_,ms_) 

diag (sin (|)m„, sin(|)m<j, O) 



(16) 



(17) 



(18) 



Making the mean-field approximation, one can obtain the 
mean-field Lagrangian as 

-Cmf = q'{-f^D^ + M}+i-f5N}~jifi)q' 

+UM+Ui$[A],$*[A],T), (19) 



where 



M} = m/+ - 4Gsf7} + 2GDi<y'f,a'f„ - if^Vf"), (20) 
N'f = mf^~AGsv'f~2Gu{o'f,v'f.+v'f'o'f„) (21) 

for/^/', /^/", f'^r and 

C/Af = 2Gs {o'/ + v'/)~^G^<t'M 

f—u,d.,s 

+4Gu{alvWs + vyWs + v'uVdO- (22) 

Performing the path integral over the quark field, one can ob- 
tain the thermodynamic potential f] (per volume) for finite T 
and /x: 



n = -2 J2 

f—u,d.s 



(2^ 



NcEf 



with Ef = a/p2 + M'f + N' 



The three-dimensional cutoff for the momentum integration 
is introduced [ 341, since this model is nonrenormalizable. For 
simplicity we assume the isospin symmetry for the u-d sec- 
tor: mi = niu = md- This three-flavor PNJL model has five 
parameters Gs, Gd, mi, nis and A. One of the typical param- 
eter sets is shown in Table These parameters are fitted to 
empirical values of pion decay constant and tt, K, 77' meson 
masses at vacuum. 



m;(MeV) 


m«(MeV) 


4(MeV) 


Gsyl' Go A" 


5.5 


140.7 


602.3 


1.835 12.36 



TABLE II: Summary of the parameter set in the NJL sector taken 
from Ref. HI]. 



For imagina ry u , Q is invariant under the extended Z3 
transformation 14411 . 



(24) 



with integer fc. This invariance ensures the Roberge-Weiss pe- 
riodicity |57] in the imaginary chemical potential region [441]. 
Any reliable model should have this extended Z3 symmetry, 
when imaginary ji is taken in the model. This is a good test 
for checking the reliability of the model. The PNJL model has 
the extended Z3 symmetry [IMIill- 

The four-quark vertex Gs is originated in a one-gluon ex- 
change between quarks and its higher-order diagrams. If the 
gluon field Ai, has a vacuum expectation value {Aq) in its time 
component, A^, is coupled to (Aq) and then to ^ through L. 
Hence we can modify Gs into an effective vertex Gs(^) de- 
pending on # ifssll . The effective vertex Gs{^) is called the 
entanglement vertex and the model with this vertex is the EP- 
NJL model. It is expected that <P dependence of Gs(^) will 
be determined in future by the accurate method such as the 
exact renormalization group method lf58l - l60ll . In this paper, 
however, we simply assume the following form for Gs (^): 



Gs{<P) = Gs[l - ai'P'P* - a2{'P^ + <P*^)]. 



(25) 



^ In [1 + 3'?e-'3(^^-^) + 3<?>*e-2^(^/-''^ + e-s/^Ci^z-P)] 
i In [1 + 3(?*e-'3(-^/+^) + 3(|>e-2/3(B/+A') _^ ^-3P(Ef+^,)^' 
Um+U{$,$*,T) (23) 



This form preserves the chiral symmetry, the charge conjuga- 
tion (G) symmetry ll48ll and the extended Z3 symmetry [44}[. 
This entanglement vertex modifies the mesonic potential Um, 
the dynamical quark masses A/j and N'^-. This is the three- 
flavor version of the EPNJL model. In principle, Gd can de- 
pend on (p, but <P dependence of Gd is found to yield qual- 
itatively the same effect on the phase diagram as that of Gs. 
Following Ref. ll35ll . we neglect <l> dependence of Gd as a 
simple setup. In the present analysis, (p dependence of Gd 
is thus renormalized in ai and a2- The EPNJL model thus 
constructed keeps the extended Z3 symmetry. 

The parameters ai and a2 in (|25] ) are fitted to two results 
of LQCD at finite T; one is the result of 2h-1 flavor LQCD at 
/i = 16 111 that the chiral transition is crossover at the phys- 
ical point and another is the result of degenerate three-flavor 
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LQCD at /X = iTw fH that the order of the RW phase transi- 
tion at the RW endpoint is first order for small and large quark 
masses and second order for intermediate quark masses. The 
parameter set (ai, 02) thus determined is located in the trian- 
gle region 

{-1.5ai + 0.3 <a2< -0.86ai + 0.32, as > 0}. (26) 

In this paper we take ai — 0.25, ^2 = 0.1 as a typical exam- 
ple, following Ref. [35]. 

The classical variables X = <P, <P*, af and 77/ are deter- 
mined by the stationary conditions 

dn/dX = 0. (27) 

The solutions to the stationary conditions do not give the 
global minimum of H necessarily. There is a possibility that 
they yield a local minimum or even a maximum. We then have 
checked that the solutions yield the global minimum when the 
solutions X{T, /i, 9) are inserted into ( |23] ). 

m. NUMERICAL RESULTS 

In this section we show numerical results for the origi- 
nal condensates (cr/, 77/,^), since this makes our discussion 
transparent. Under the parity transformation, cr/, 77/ and <P are 
transformed into , —rjf and ^, respectively. This means that 
rjf is 9 odd while aj and <P aie 9 even, since the Lagrangian is 
invariant under the combination of the parity transformation 
and the transformation 6 —6. Thus 7;/ is an order param- 
eter of the spontaneous parity breaking, while a / and <P are 
approximate order parameters of the chiral and the deconfine- 
ment transitions, respectively. As an approximate order pa- 
rameter of the chiral transition, ai = a-a = (Jd is more proper 
than (Ts, since 777/ -C rris. 

A. Thermodynamics at /i = 

In this subsection, we consider the case of /i = where 
charge conjugation symmetry (C) is exact. Meanwhile, parity 
symmetry (P) is exact only at = 0, ± tt, since e'^ agrees 
with e"''' in © when 6* = 0, ± tt. 

Figure[T]shows T dependence of cr/, as and ^ at ^? = /i = 0, 
where cr/ and (Js are normalized by uo = cr/ at T = /i = 
9 = Q. Here cr; and ^ describe the chiral and deconfine- 
ment transitions, respectively. In the PNJL model of panel 
(a), the chiral restoration transition takes place after the de- 
confinement transition. In the EPNJL model of panel (b), 
meanwhile, both the transitions occur simultaneously. In the 
EPNJL model, as decreases rapidly near the pseudocritical 
temperature Tc = 160 MeV, but goes down gradually above 
Tc. The rapid change of as comes from that of For both the 
PNJL and EPNJL models, 77/ and 77^ are zero at any T. The P 
symmetry is thus always preserved when 6' = 0. 

Figure|2]shows 9 dependence of fi and the order parameters 
at r = 7i = in the EPNJL model; note that the EPNJL 
model agrees with the PNJL model at T = 0, since Gs = 
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Fig. 1: T dependence of the order parameters at S = and ^ = Q 
in (a) the PNJL model and (b) the EPNJL model. The solid, dashed 
and dotted lines represent cr;, Oa and respectively, where cr; and 
Os are normalized by ao — ai{T = jj, — 6 = 0) . 



Gs there because of<P = 0. As shown in panel (a), f2is9 even 
and has a cusp at 6* = tt. This indicates that a first-order phase 
transition takes place at T = f.i = and 6* = vr. As shown 
in panel (b), meanwhile, the r]f aie 9 odd, while a/ and <P 
are 9 even. The condensate 77/ and 77.; have jumps at 9 = tt, 
indicating that the first-order transition mentioned above is the 
spontaneous parity breaking. This is nothing but the Dashen 
phenomena [18]. 

Figure [3] shows 9 dependence of the order parameters and 
the effective quark mass 77/ ee y^A/j.^ + N'^.'^ at T = 

163 MeV and = in the EPNJL model. For this higher tem- 
perature, the Dashen phenomena do not take place at 6* = vr. 
Actually r/i and r]s vanish there, although they become finite 
at9 0, TT, 2tt where P is not an exact symmetry. The other 
order parameters, aj and are smooth periodic functions of 
9. The Polyakov loop becomes maximum at 6* = tt, since 
the effective quark mass 77/ becomes minimum and the ther- 
mal factor exp{—l3Ef) is maximized in (l23T l. 

Figure |4] shows T dependence of the order parameters at 
9 ^ TT and fi ~ 0. Comparing this figure with Fig.[T] one can 
also see 9 dependence of the order parameters. In the PNJL 
model of panel (a), |77i| and It;^! are finite below the critical 
temperature Tp — 194 MeV and vanish above Tp. Thus the 
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Fig. 2: 9 dependence of (a) fl and (b) the order parameters at T = 
H = Q m the EPNJL model. In panel (a), J7o = fl{0 = 0) is 
subtracted from O. See the legend for the definition of lines. 
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-0.015 




Fig. 3: 6 dependence of (a) the order parameters and (b) the effective 
quark mass 11 f at T = 163 MeV and ^ = in the EPNJL model. 
In panel (b), 11 f is normalized by the value atr = /i = 6 = and 
the normalized 77/ is compared with the Polyakov loop See the 
legend for the definition of lines. 



P symmetry is broken at smaller T, but restored at higher T. 
In the two-flavor PNJL model, this P restoration is second 
order llsoll . This is the case also for the present 2+1 flavor 
PNJL model. The second order P restoration induces cusps 
in \(Ti \ and \as\ when T = Tp, although the cusp is weak in 
CTs |. This propagation of the cusp can be understood by the 
extended discontinuity theorem of Ref. ll47ll . In the EPNJL 
model of panel (b), the P restoration occurs at Tp — 158 MeV 
as the first-order transition. The same property is seen in the 
two-flavor EPNJL model The first-order P restoration 
generates gaps in \ai\ and \as\ when T = Tp, although the 
gap is tiny in \as\. This propagation of the gap can be under- 
stood by the discontinuity theorem by Barducci, Casalbuoni, 
Pettini and Gatto 1,63,1 . Thus the Dashen phenomena are seen 
only at lower T, and the order of the P violation at the critical 
temperature Tp depends on the effective model taken. 

Theoretical prediction on the critical temperature of the chi- 
ral transition at 6* = and fi = and the P restoration atd = tt 
and /i = is tabulated in Table Hill Atd = 0, the chiral transi- 
tion is crossover in all of the NJL, PNJL, and EPNJL models 
At 9 ~ TT, the order of the P restoration is first order in the 
EPNJL model, but it is second order in the PNJL and NJL 
models. 



Model 9 = 9 = n 



NJL 


177 (crossover) 


170 (2nd order) 


PNJL 


200 (crossover) 


194 (2nd order) 


EPNJL 


162 (crossover) 


158 (1st order) 



TABLE III: Theoretical prediction on the critical temperature of the 
chiral transition at 9 = and fi = and the P restoration at 6* = tt 
and /I = 0. The values are shown in units of MeV. 



B. Thermodynamics at /i > 

In this subsection, we consider the case of /j, > where 
C symmetry is not exact. In general, the relation <P = <P* is 
not satisfied for finite /i, although <P and <P* are real fsO]. This 
situation makes numerical calculations quite time-consuniing. 
However, the deviation (p — tp* is known to be very small 115011 . 
For this reason, the assumption <p = <p* has been used in 
many calculations. Therefore we use the assumption also in 
this paper. 

Figure |5] represents T dependence of the order parameters 
at 6* = TT and /i = 300 MeV in the PNJL and EPNJL models. 
The P restoration takes place at high T, since rji and rj^ are 
zero there. The critical temperature of the P restoration is 
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Fig. 4: T dependence of the order parameters at 6 = tt and /i = 
in (a) the PNJL model and (b) the EPNJL model. See the legend for 
the definition of lines. 



Fig. 5: T dependence of the order parameters at ^ = tt and fj, = 
300 MeV in (a) the PNJL model and (b) the EPNJL model. See the 
legend for the definition of lines. 



Tp = 110 MeV for the PNJL model and Tp = 99 MeV 
for the EPNJL modeL For = 300 MeV, the order of the 
P restoration at T = Tp is first order in both the PNJL and 
EPNJL models. Thus the quark-number chemical potential /i 
lowers Tp and makes the P restoration sharper. 

Figure |6] shows the phase diagram of the chiral transition 
in the fjL-O-T space. The diagram is miiTor symmetric with 
respect to the fi-T plane al — Q, so the diagram is plotted 
only at 6* > 0. Panels (a) and (b) correspond to results of the 
PNJL and EPNJL models, respectively. In the ^-T plane at 
{) < 6 < TT, the solid line stands for the first-order chiral tran- 
sition, while the dashed line represents the chiral crossover. 
The meeting point between the solid and dashed lines is a crit- 
ical endpoint (CEP) of second order. Point C is a CEP in the 
^-T plane at 61 = lH H . For both the PNJL and EPNJL 
models, the location of CEP in the ^-T plane moves to higher 
T and lower jj, as 9 increases from to tt. 

In the /i-T plane at 6 = n, P symmetry is exact and hence 
we can consider the spontaneous breaking of P symmetry in 
addition to the chiral transition. For the PNJL model of panel 
(a), both the first-order chiral transition and the first-order P 
restoration take place simultaneously, and the second-order P 
restoration and the chiral crossover coincide with each other. 
The first-order and the second-order P transition line are de- 
picted by the solid and dashed lines, respectively. The meeting 



point A is a tricritical point (TCP) of the P-restoration transi- 
tion. For the EPNJL model of panel (b), the chiral and the P 
restoration transition are always first order and hence there is 
no TCP 

In the PNJL model of panel (a), the dotted line from point 
C to point A is a trajectory of CEP as 9 increases from to 
TT. Thus the second-order chiral transition line ends up with 
point A. This means that the CEP (point C) at 6* ~ is a 
remnant of the TCP (point A) of P restoration at 6* = tt. In the 
EPNJL model of panel (b), no TCP and then no CEP appears 
in the fi-T plane at 6* = tt. The second-order chiral-transition 
line (dashed line) starting from point C never reaches the fi-T 
plane at 61 = tt. 

Figure |7] snows the projection of the second-order chiral- 
transition line in the fi-9-T space on the fi-9 plane. The 
solid (dashed) line stands for the projected line in the EP- 
NJL (PNJL) model. The first-order transition region exists on 
the right-hand side of the line, while the left-hand side coiTe- 
sponds to the chiral crossover region. The first-order transition 
region is much wider in the EPNJL model than in the PNJL 
model. In the EPNJL model, eventually, the chiral transition 
becomes first order even at ^ — when 9 is large. 
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Fig. 6: Phase diagram of the chiral transition in the fi-6-T space. 
Panel (a) shows a result of the PNJL model and panel (b) corresponds 
to a result of the EPNJL model. 
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Fig. 7: The projection of the second-order chiral-transition line in 
the IJ.-9-T space on the ji-O plane. See the legend for the definition 
of line. 



C. The sign problem on LQCD with finite 6 

In the PNJL Lagrangian ( fTSI l after the transformation (|2]), 
6 dependence appears only at the light quark mass terms, 
771/ cos(6'/2) and 7T7; sin(6'/2). These terms are much smaller 
than Aqcd as a typical scale of QCD. This means that the 
condensates, ctj, a'^, r/l and 77I,, have weak 6 dependence. This 
statement is supported by the results of the PNJL calculations 
shown in Fig. [8] 

The sign problem is induced by the 9 odd 777; sm{9 /2) term. 
The 6'-odd (P-odd) condensates, rj^ and 77^, are generated by 
the ^-odd mass term. One can see in Fig. [8] that the 6'-odd 
condensates are much smaller than the 6-even condensates, ct; 
and a'g. This fact indicates that effects of the 6'-odd mass term 
are negligible. Actually, if the term is neglected, the 6'-even 




Fig. 8: 6 dependence of the order parameters, ctj, a'^, rjl and rj'g, at 
T — fi — calculated with the EPNJL model. See the legend for the 
definition of lines. 



condensates change only within the thickness of line, while 
the ^-odd condensates vanish. The neglect of the ^-odd mass 
is thus a good approximation. 

The validity of the approximation can be shown more ex- 
plicitly in the following way. The 6'-odd (P-odd) condensates, 
rj'i and 77^., are zero at 6* = 0, since the 6'-odd mass vanishes 
there. The weak 9 dependence of 77^ and 77^ guarantees that rj'i 
and are small for any 9. Setting 77; =1/^=0 in All '^^'^ 
leads to 



A// = 



cos 



sm 



;i)77i,, 



(28) 
(29) 



where Af/ « /Iqcd and N'l « mi. Since the thermodynamic 

the term A^/^ is negligible 



/2 



potential is a function of 
compared with Af/^. 

In LQCD, the vacuum expectation value of operator O is 
obtained by 

(O) = j VAO{detMi{9)f dctMse-^o (30) 
= J VAO' {dctM'i{9)f detMsC-^' (31) 
with the gluon part Sg of the QCD action and 



^, ^ ^{dctMi{9)y 



idetM[{9)y 



(32) 



where detM[{9) is the Fermion determinant in which the 
9-odd mass is neglected and hence has no sign problem. As 
mentioned above, one can assume that 



detMi{9) 
detM'i{9) 



1. 



(33) 



Thus the reweighting method defined by (ISTT l may work well. 
In the 0-even mass, 771; cos(0/2), the limit of 9 ~ tt corre- 
sponds to the limit of 777; = with 777^ fixed. Although the 
limit is hard to reach, one can analyze the dynamics at least at 
small and intermediate 9. 
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IV. SUMMARY 

We have investigated effects of the theta vacuum on the 
QCD phase diagram for the realistic 2+1 flavor system, us- 
ing the three-flavor PNJL and EPNJL models. The effects can 
be easily understood by the S'C/a(3) (S) C^a(1) transformation 
(|2]l. After the transformation, the 6'-odd mass, m; sm{9/2), 
little affects the dynamics, so that the dynamics is mainly gov- 
erned by the 0-even mass, mi cos(6'/2). In the 9-even mass, 
the increase of 9 corresponds to the decrease of m; with 
fixed. This means that the chiral transition becomes strong 
as 9 increases. This is true in the results of both PNJL and 
EPNJL calculations. Particularly in the EPNJL model that is 
more reliable than the PNJL model, the transition becomes 
first-order even at ^ = when 9 is large. This result is im- 
portant. If the chiral transition becomes first order at /i = 0, 
it will change the scenario of cosmological evolution. For ex- 
ample, the first-order transition allows us to think the inhomo- 



geneous Big-Bang nucleosynthesis model or a new scenario 
of baryogenesis. 

Using the fact that the 6'-odd mass is negligible, we have 
proposed a way of circumventing the sign problem on LQCD 
with finite 9. The reweighting method by defined dSTT i may 
allow us to do LQCD calculations and get definite results on 
the dynamics of the 9 vacuum. 
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